Intermediate range chemical ordering of cations in simple molten alkali halides 
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The presence of first sharp diffraction peaks in the partial structure factors is investigated in 
computer simulations of molten mixtures of alkali halides. An intermediate range ordering appears 
for the Li + ions only, which is associated with clustering of this species and is not reflected in the 
arrangement of other ions. This ordering is surprising in view of the simplicity of the interionic 
interations in alkali halides. The clustering reflects an incomplete mixing of the various species on a 
local length scale, which can be demonstrated by studying the complementary sub-space of cations 
in the corresponding pure alkali halides by means of a void analysis. 



Prepeaks or first sharp diffraction peaks in the diffrac- 
tion structure factors of amorphous solids or liquids are 
associated with intermediate-range order- 1 . In simple liq- 
uids, like hard-spheres or simple MX molten salts the 
structure factors S(k) are dominated by a principal peak 
at k pp ~ 2zc, where a is the nearest neighbour separa- 
tion, and have only a small amplitude for smaller values 
of k, monotonically decaying as k — > 0-. This effect, 
which is captured in simple integral theories of such flu- 
ids, means that, although the radial distribution func- 
tions of these systems do exhibit oscillatory structure 
out to much larger separations than a, it is simply a 
consequence of the packing of the particles reflected in 
the nearest-neighbour separation. Alternatively, viewed 
in reciprocal space, the density fluctuations reflected in 
S(k) are suppressed for k < k pp by packing effects in 
van der Waals atomic liquids and by the combination of 
packing and coulomb ordering in simple molten salts. 

Exceptions are found when a competing interaction 
affects the order. In the atomic system phosphorus^, 
where the competition between tendency to form lo- 
cal quasi-molecular tetrahedra, which determines the 
nearest-neighbour separation, and the packing of these 
tetrahedra in a disordered network introduces new struc- 
ture on the intermediate lengthscale and results in a 
substantial prepeak. Prepeaks have been most exten- 
sively studied in the partial structure factors of network- 
forming MX2 binary systems, like SiO^ and ZnCh 5 
where, again, they reflect a competition between the for- 
mation of local coordination polyhedra (in these cases, 
tetrahedral SiC>4 and ZnCU units) and the arrangement 
of these units in three dimensions which is influenced 
by the preferred value for the cation-anion-cation bond 
angle 6,7 . Prominent prepeaks are also seen in MX3 sys- 
tems, for similar reasons^. In these cases it has been 
found useful to associate the position of the prepeak with 
the arrangement of voids in the spatial distribution of the 
species of interest^; the ordering of the voids reflects 
the inhomogeneities introduced in the arrangements of 
atoms due to these competing structural influences. The 
void-void structure factor shows a principal peak which 
coincides with the position of the prepeak in the atomic 
structure factor. In the simple systems which lack pre- 



peaks, the void and atomic principal peaks coincide. 

In this article we are concerned with the extent to 
which such ideas can be transferred to discuss the struc- 
ture of ternary mixtures (and also of binary mixtures of 
metals which might also be thought of as ternary, with 
the valence electrons thought of as the third structural 
component). Here "prepeaks" may be observed for essen- 
tially trivial reasons, and such cases should be discounted 
from consideration. If, for example, the network-former 
ZnCl2 is mixed with a network-breaker, like RbCl, to 
form a ternary system, the network is broken up into 
ZnCl 2 ~ molecular ion units if more than two moles of 
RbCl are combined with one of ZnCki. A first sharp 
diffraction peak is now seen in the Zn-Zn partial struc- 
ture factor— but it arises because the liquid should be 
considered as a simple unassociated mixture of Rb + ions 
and ZnCl 2- ions and the "prepeak" can be associated 
with the principal peak of the partial structure factor of 
the latter—. This structure can be understood from sim- 
ple molten salt theories. On the other hand, in appar- 
ently similar cases, such as when the structure-breaker 
Na2 is added to SiC>2 at equal mole fractions (thereby 
allowing only incomplete network breakdown) the Na- 
Na partial structure factor shows a prepeak and the 
associated lengthscale corresponds to the distance be- 
tween Na-containing channels in the remnants of the 
SiC>2 network-^. Since these channels are associated with 
Na conduction in the glassy state, in this case the struc- 
ture evident in the prepeak is clearly associated with very 
significant effects on the ion dynamics. More general con- 
siderations link non-trivial relationships between the sin- 
gle particle and collective dynamics to the existence of 
a prepeak in the structure factor—, and these consider- 
ations motivate, in part, our interest in the alkali halide 
systems. 

In the cases mentioned above, the appearance of (non- 
trivial) prepeaks is associated with strong interactions 
involving at least some components. However, there is 
one system in which a prepeak has been reported where 
one would have thought all the interactions were quite 
simple. Ribeiro^ reported a prepeak in the partial Li-Li 
structure factor in the alkali halide mixtures LiF:KF and 
LiChKCl from computer simulations. Here both com- 
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FIG. 1: F~-F~, K+-K+ and Li+-Li+ partial structure factors 
for various compositions of LiF:KF mixtures. Sff(A0 and 
<5 , KK(fc) were respectively shifted by values of 3 and 1.5. 



ponents are MX molten salts, whose structures are well 
understood on the basis of packing and coulomb order- 
ing, as recaptured in simple HNC theories for example 2 . 
We will examine these and several related alkali halide 
mixtures in the present work. Our objective is to un- 
derstand how the intermediate-range order arises in this 
system and to relate it to that seen in other, more famil- 
iar examples. 

In his simulations, Ribeiro used rigid ion interaction 
models for all the ions. In recent years we have shown 
that taking polarization effects into account is impor- 
tant to obtain a description of the physical properties 
of molten salts-^ which is transferable between the pure 
salts and their mixtures. Even for the alkali halides, 
where effective pair potentials reproduce many proper- 
tics remarkably well, polarizable potentials are necessary 
to reproduce ab initio derived forces accurately^ and 
therefore to faithfully represent the real interactions. We 
therefore undertook new simulations of mixtures of alkali 
halides using polarizable ion models which have been de- 
rived by force-fitting from first-principles—, extensively 
tested, and proven to give a good representation of the 
materials. 

We have performed MD simulations of KF, NaF and 
KC1 pure melts and of LiF:KF {x KF =0.25, 0.50, 0.75) 
LiF:NaF (x NaF =0.75) and LiChKCl (x KC i =0.41) mix- 
tures in the canonical ensemble. The simulation cells 
contained a total of 432 ions, or 3456 ions when neces- 
sary, and were equilibrated in the NPT ensemble at a 
zero pressure. The corresponding densities were in close 
agreement with the experimental ones. 

The like-like partial structure factors are given on fig- 



ure [T] for the several compositions of LiF:KF mixtures. 
For K + and F _ ions, only one very well defined (princi- 
pal) peak is observed and it occurs close to k pp ~ — 
where a is the nearest-neigbour separation identified 
from the corresponding partial radial distribution func- 
tion. On the other hand, S-Li-u(k) displays two peaks for 
all the compositions. The prepeak (or FSDP) is a broad 
peak and its intensity rises significantly with the KF con- 
centration of the mixture. Even if the broadness does not 
allow us to determine an exact position for that peak, it 
clearly shifts significantly with the change in composi- 
tion, passing from kpsDP ~ 0.82 A -1 for xkf = 0.25 
to kpsDP ~ 1-18 A -1 for ikf = 0.75. Several temper- 
atures ranging from 773 K to 1300 K were investigated. 
The corresponding number density of ionic pairs shifts 
from 2.94 x 10~ 2 A~ 3 (x KF = 0.50 and T=773 K) to 
2.04 x 10" 2 A" 3 (i K F = 0.75 and T=1300 K). For a 
given composition, the positions and widths of the FS- 
DPs remained unchanged from one temperature to an- 
other. The peak at higher k corresponds to the principal 
peak and occurs at k pp ~ — . The lithium ions spatial 
distribution is thus characterized by the presence of a sec- 
ond length scale which does not seem to affect the other 
ions. These findings closely parallel Ribeiro'si 3 -. 

This behavior indicates the existence of chemical or- 
dering. In a binary system, topological versus chemi- 
cal ordering can be investigated by using the Bhatia- 
Thornton formalis m 16 ' 17 . Three partial structure fac- 
tors are constructed which correspond to the number- 
number (5jvjv(fc))i number-concentration (Sjvc(&)) and 
concentration-concentration (Scc(k)) fluctuations. In 
a ternary mixture, the same definition can be used for 
S NN (k): 



(k) 



a,/3 



Cf3S a p(k), 



(1) 



where c a is the ionic fraction of chemical species a. 
As regards Scc(k), a general expression for a multi- 
component fluid was developed by Bhatiai^, but from 
our perspective it is more chemically intuitive to study 
fluctuations for a given pair of atoms in the same way as 
in a two component fluid. We have therefore determined 
the following partial Sc a c {k) structure factors: 

Sc a C (k) = c a cp(S aa (k) + Sp^k) - 2S a p(k)). (2) 

The two concentration-concentration functions involving 
Li + ions and the number-number one have been plotted 
on figure [5] They give us a first insight into the ori- 
gin of the intermediate-range ordering observed for that 
species. The only function which presents a peak in the 
low-A; region is 5 , c Li c K (^)- This means that the Li + ions 
are chemically ordered with respect to the K + but not to 
the F - ions. 

Unlike the tetrahedral network forming ionic liquids 
like Si02 or ZnCki, no well-defined first coordination shell 
structure is formed around the lithium cations and their 
fluoride ion coordination number can be either 3, 4 or 5. 
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FIG. 2: Bhatia-Thornton structure factors for the LiF:KF 
mixture at ikf = 0.50 composition. 



Therefore, the examination of the influence of a partic- 
ular ion solvation shell on the structures occurring at a 
medium range scale does not give any clue to the origin 
of the intermediate range ordering. On the contrary, vi- 
sual examination of the Li + ion positions in individual 
configurations of the liquid reveals an heterogeneous dis- 
tribution as shown on figure [3] Regions of high and low 
concentrations can be identified, showing that the second 
length scale observed in Su-Li(k) is associated with the 
formation of clusters of Li + ions. 

Such a clustering means that the local concentration in 
Li + ions fluctuates in the liquid, revealing an incomplete 
mixing of LiF and KF at the nanoscopic scale. As a con- 
sequence, the K + local concentration might be expected 
to fluctuate too, but the absence of low-fc features in the 
like- like partial structure factor of that ion shows that its 
distribution is homogeneous. What emerges is a picture 
of a homogeneous but disordered (i.e. non network- like) 
matrix consisting of K + and F~ ions in which clusters of 
Li + ions are formed. Note that the position of the prin- 
cipal peak in Skk is not shifting. The behaviour seen 
resembles closely that found by Voigtmann et al in re- 
cent experimental studies of the binary metallic liquid 
Ni:Zr— . In that case, the Ni-Ni partial structure showed 
a double peak structure of very similar character to that 
of the Li ions in our calculation, whereas the larger host 
Zr ions behaved analogously to the K + ions in our cal- 
culations: they showed a normal simple metallic struc- 
ture factor which was qualitatively reproduced by that 
of hard-spheres. 

In simulations of pure KF, the K + ions can be arbitrar- 
ily assigned to two groups, which are labelled as "hosts" 
and "guests" , the relative proportions of the two groups 
can be varied in order to correspond to the compositions 
of the LiF:KF mixtures. It is then of interest to study the 
possible positions at which another species could be in- 
corporated into the melt if the K + guests were removed. 
To do this, it is possible to examine the spatial arrange- 
ment of the complementary space of the hosts ions, i.e. 




FIG. 3: Snapshot of the Li ion positions in a LiF:KF mixture 
(ikf = 0.75). The whole cell is projected onto a plane. Pairs 
of Li + ions separated by a distance smaller than the Li + -Li + 
RDF first minima were connected by "bonds". Note that 
because of the boundary conditions, some of these bonds are 
not appearing here. 



of the regions where those ions are absent. Such "voids" 
are not physical voids as they include F~ anions and K + 
guests. They can be identified by performing a Voronoi 
analysis of the host cation positions in the instantaneous 
liquid configuration—. The vertices of the Voronoi poly- 
hedra in a disordered structure define the center of a 
group of four atoms (here four host cations), which are 
mutually nearest neighbors, and are known as the Delau- 
nay simplices (DS) of the structure. Around each DS a 
circumsphere which passes through the four atoms may 
be drawn and, since no other atom center lies within the 
sphere, the radius r c of the circumsphere gives a measure 
of the empty (void) space in between the atoms. From 
the positions of the center of each sphere, a void- void 
structure factor can be calculated: 

Svv(fc) = E ex P( ik ' Rii )^ • ( 3 ) 

Note that Syv(k) does not correspond to the partial 
structure factor of the "guest" K + ions in the pure melt 
simulation. Since these K + guests were assigned arbi- 
trarily (i.e. randomly), their structure factor is simply 
proportional to SkkQz) for the pure melt (plus an inco- 
herent background). 

In order to study the complementary space of the host 
K + ions in pure KF, the positions of the voids between 
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FIG. 4: Void-void (top curves) and Li + -Li + (bottom curves) 
partial structure factors in LiF:KF. The void-void structure 
factors were determined in pure KF, where the K + ions have 
been separated into "hosts" and "guests" . The voids positions 
were determined from the positions of the "host" ions only 
(see the text). The left hand ordinate axis corresponds to the 
void-void structure factors and the right hand ordinate axis 
to the Li + -Li + ones. 



them have been extracted from our trajectories, and the 
corresponding structure factor was computed. The re- 
sults obtained for several proportions of hosts K + ions 
are shown on figured! where the Li + -Li + partial struc- 
ture factors of the corresponding LiF:KF compositions 
are also given. One can see that the positions of the main 
peak of Syy(k) and of the FSDP of Su-u(k) coincide 
for all the compositions. This shows that the interme- 
diate range ordering of the Li + ions in LiF:KF occurs 
at the characteristic length scale of the complementary 
sub-space of the corresponding number of K + ions in pure 
KF, which is consistent with the view of LiF:KF mixtures 
as a partial homogeneous KF matrix of which voids are 
filled with Li + ions. 

It is natural to ask if this intermediate range struc- 
turing is also present in other simple mixtures of al- 
kali halides like LiF:NaF and LiChKCl. The partial 
structure factors computed for those systems showed ex- 
actly the same features as in LiF:KF mixtures, i.e. the 
presence of a single peak in the anion structure fac- 
tors and in the structure factors involving both K + and 
Na + ions, whereas Su-u(k), which is shown on fig- 
ure presents two peaks. In LiChKCl, the FSDP is 
at kpsDP ~ 0.98 A -1 and in LiF:NaF it appears at 
kFSDP ~ 1-31 A -1 . Once again, the positions of those 
FSDP coincide exactly with the principal peak of the void 
void structure factor of the corresponding pure host melt 



(KC1 or NaF). This shows that the same origin for inter- 
mediate range ordering of Li + ions is observed in all the 
mixtures and that its associated length scale is dependent 
on the other component of the mixture. The comparison 
of LiChKCl and LiF:KF mixture also edves iniDortant in- 
i t- 
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FIG. 5: Void-void (top curves) and Li + -Li + (bottom curves) 
partial structure factors. The void-void structure factors were 
determined in pure KC1 and NaF, where the K + or Na + ions 
have been separated into "hosts" and "guests". The voids 
positions were determined from the positions of the "host" 
ions only (see the text). The SLi-Li(fc) were computed in 
LiCbKCl and LiF:NaF mixtures. The left hand ordinate axis 
corresponds to the void-void structure factors and the right 
hand ordinate axis to the Li + -Li + ones. 



those melts, and so is the Li + -Li + first neighbor distance 
(given by the principal peak of Su-u{k)) ■ On the con- 
trary the intermediate range ordering of Li + ions appears 
at similar length scales, which means that the knowledge 
of the first shell structure does not help in understanding 
the clustering properties of the melt. 

The intermediate range chemical ordering of Li + ions 
then is a general property of simple mixtures of lithium 
halides with the other alkali halides. It is due to poor 
mixing properties on the nanoscopic scale, and the liq- 
uid formed consists in a matrix of "host" ions filled with 
clusters of Li + . The length scale associated with those 
clusters is set by the void- void correlation in the pure host 
liquid, and no indication is given by the real-space short- 
range structure. This result shows that the formation of 
a network may not be a prerequisite to the formation of 
those clusters, and hence of intermediate range chemical 
ordering in condensed matter systems. 
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